Itinerant vs Localized Heavy-Electron Magnetism 
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It is demonstrated that itinerant-localized transition of heavy electrons occurs inside the magnet- 
ically ordered phase of the Kondo-Heisenberg lattice. The phase diagram and electronic structure 
are derived by means of the continuous-time quantum Monte Carlo combined with the dynamical 
mean-field theory. Around the itinerant-localized transition, nearly flat bands appear on the Fermi 
surface with almost vanishing quasi-particle renormalization factor. At the same time, there emerges 
a strong local magnetic fluctuation with a minute energy scale. Considering both antiferromagnetic 
and ferromagnetic Heisenberg interactions, coherent understanding is achieved on rich phase dia- 
grams observed in CeRhi-a; Covins, CeRu2(Si a: Gei_ a ;)2, UGe2 and CeT2Alio (T=Fe,Ru,Os). 



It has been a long-standing problem whether and how 
nearly localized / electrons change their character de- 
pending on the environment. The localized / electrons 
tend to form magnetically ordered states with intersite 
exchange interactions, while they acquire the itinerancy 
by the Kondo effect resulting in a heavy Fermi liquid. 
Intriguing phenomena such as quantum criticality and 
unconventional ordered states have been observed in the 
competing region, and hence understanding the distinc- 
tion or duality between itinerant and localized characters 
is a fundamental and important issue in strongly corre- 
lated electron systems. 

The difference between itinerant and localized charac- 
ters appears in the Fermi surface; itinerant / electrons 
contribute to the Fermi volume, while localized ones do 
not. It has been shown for the Kondo lattice in large 
dimensions [14| that / electrons remain itinerant at the 
magnetic quantum critical point (QCP). In this paper, 
therefore, we focus on possible change of /-electron char- 
acters inside magnetically ordered states. Experimen- 
tally, CcRhi_ :r Co :r In5 |l| shows, for example, that the 
Fermi surface at x = resembles the case for localized 
/ electrons, while at x ~ 0.4 which is deep inside the 
antiferromagnetic (AFM) phase [2|-|4|, the Fermi surface 
corresponds to the itinerant one. The itinerant-localized 
transition (ILT) in the ordered state has also been found 
in CeRu2(Si a; Gei_ a ;)2 accompanying the magnetic tran- 
sition between ferromagnetic (FM) and AFM states |5J. 
Furthermore, drastic Fermi-surface reconstruction inside 
the FM phase has been observed in a heavy-electron su- 
perconductor UGe2 [a-lSl- Thus it is desirable to explain 
these behaviors in a unified manner, which should lead 
to deeper understanding of the itinerant-localized duality 
of electrons. Obviously it is necessary to go beyond the 
celebrated Doniach-type phase diagram. 

Many theoretical attempts have already been made to 
explain the ILTs [9h14J . The presence of ILTs inside the 
ordered phase has been pointed out within Gutzwiller- 
type variational [15| . and mean- field |16j approximations. 
These approximations, although simple, are unable to 
deal with the duality of electrons; the itinerant and local- 
ized behaviors depend on the energy scale. On the other 



hand, numerical approaches must be accurate enough to 
deal with a minute energy scale that can be even smaller 
than the Kondo temperature in the quantum critical re- 
gion. In this pap er, we use the continuous-time quantum 
Monte Carlo [121 llSl combined with the dynamical mean- 
field theory (DMFT) [ljj H3| which takes full account of 
local correlations, and becomes exact in the limit of high 
dimensions. As we shall show in the following, the ILT 
of / electrons is controlled by sign and magnitude of the 
Heisenberg exchange. The critical temperature can be 
tuned to minute, which makes the first-order transition 
to a quantum critical transition. There appears consid- 
erable asymmetry on both sides of the ILT that can be 
probed by a large effective mass and enhanced local mag- 
netic susceptibility. The resultant phase diagram leads to 
a coherent understanding for the experimentally observed 
ILTs explained above. 

In order to discuss the ILT, we consider the following 
Kondo-Heisenberg lattice: 
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where £& = £& — \x with /x being the chemical poten- 
tial. The first and second terms are the kinetic energy of 
conduction (c) electrons and local Kondo interaction, re- 
spectively. The Heisenberg interaction in the third term 
realizes the simplest localized magnetism, and is incor- 
porated in the DMFT as the mean field. The summa- 
tion with respect to (if) is taken over pairs of the near- 
est neighbor sites with the coordination number z. We 
take a bipartite lattice with the semi-circular density of 
states D(e) = {A/ttW 2 )^W 2 - 4e 2 . The band width 
W = 1 is the unif of energy. Even if the inter-site in- 
teraction Jh is small as compared to Jk, Jh can still 
be important since it competes with the Kondo energy 
7k oc exp[— l/D(/j,) Jk] instead of Jk, as will be demon- 
strated later. Throughout this paper, the number of c 
electrons is fixed as n = 0.95 per site. 

Let us first derive the characteristic energy scale from 
temperature dependence of the electrical resistivity p(T). 
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FIG. 1: (Color online) Temperature dependence of the elec- 
trical resistivity for (a) Jk > Jci m 0.16 and (b) Jk < Jc2, 
where J C 2 gives the ILT as discussed in the text. The dotted 
arrow shows the change of Tp with decreasing Jk. 



Wc use the formula 



P(TY 



a lira ImII(q = Q,u>)/ui. 



(2) 



Here n(q, u>) is the current-current correlation function 
and a is a constant. Since the local vertex correction in 
the DMFT does not contribute to the conductivity |20j], 
we evaluate a simple particle-hole bubble. 

Figure [1] shows the temperature dependence of resis- 
tivity for the various values of Jk- We fix the Heisenberg 
exchange as Jh = 0.025 except for Fig.|U For Jk = 0.350 
shown in Fig. QJa), the system remains paramagnetic 
down to zero temperature. The resistivity at high T in- 
creases with decreasing T due to the impurity-like Kondo 
effect. On the other hand, the resistivity shows the 
metallic behavior at low temperatures reflecting the de- 
velopment of coherence. The characteristic temperature 
Tp ~ 0.02, which gives the peak of p(T), corresponds to 
the effective Fermi temperature below which the heavy 
Fermi liquid is formed. As seen in Fig. QJa) , Tp decreases 
with decreasing Jk- 

With J K < Jci ^ 0.27, the system undergoes an AFM 
transition at the Neel temperature Tn as indicated by ar- 
rows in Fig. [TJ The resistivity shows almost no anomaly 
at the transition point with Jk = 0.250. The kink at 
Tn becomes clearer for smaller Jk- Contrary to the case 
with larger Jk shown in the upper panel (a), the effec- 
tive Fermi temperature increases with decreasing Jk for 
smaller Kondo interactions as shown in Fig. [TJb). 

We plot in Fig. [2ja) the Neel temperature Tn and 
Fermi temperature Tp as functions of Jk- Note that Tp 
is finite at the magnetic QCP with Jk = J c i — 0.27, but 
it becomes almost zero at Jk = J C 2 — 0.16 inside the 
AFM phase. This point actually separates the itinerant 
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FIG. 2: (Color online) (a) Temperature (T) vs Kondo inter- 
action (Jk) phase diagram for Jh = 0.025. The Jk depen- 
dences of (b) magnetic moment, (c) renormalization factors 
and (d) local magnetic susceptibility are also shown. 



and localized behaviors of / electrons, and is called the 
ILT point. More details will be discussed later around 
Fig. [3] From the localized side, Tp tends to zero toward 
J C 2 much more rapidly as compared to the itinerant side. 

Figures [2jb), (c), (d) show relevant physical quanti- 
ties around the ILT. The spontaneous magnetic moment 
M = |(5? + s* 4 )| in (b) shows a kink at the ILT point 
at T = 0.001, while it becomes smooth at T = 0.005. 
This behavior reflects the first order transition at tem- 
peratures lower than T = 0.001. Namely the critical 
end point where the first-order transition line terminates 
should be located at certain T < 0.001. We have not ob- 
served the ILT at half filling with n = 1 where the system 
is insulating. Hence, the temperature at the critical end 
point should tend to zero asn^ 1. 

Figure EJc) shows the renormalization factors a\f, 
which is defined by 



as ACT 
^Tol 1 - — 



(3) 



for sublattice components A = A,B and spin a =t,4- 
Here E>, cr (e) is the local self energy of c electrons. Note 



that the relations £a-|- = ^Bi and Eaj, = ^Bt hold 
in the AFM phase. We evaluate a\ a numerically from 
the low-energy behavior of the self energy at sufficiently 
low temperatures. For example, we take T = 0.001 for 
J K = 0.225 and T = 0.0025 for J K = 0.35. In the param- 
agnetic state, oa-j- (= aet) decreases with decreasing Jk 
reflecting the reduction of 7k. 

In the AFM phase, ast is strongly reduced near the 
ILT point, while OAf is enhanced toward the ILT from the 
itinerant side. Note that we have the relation riAt -> 71 Bt 
with n\f being the number of c electrons with f spin at 
sublatticc A. In the localized side, oa^ i s almost unity. 
Namely, the large renormalization takes place only for c 
electrons with the minority spin at each site. As shown 
elsewhere [2l| , the geometric mean y/aXfOBt contributes 
to the specific heat coefficient for simple bipartite lattices. 
This quantity is also plotted in Fig. [2tc) , and behaves in 
a manner similar to aef- The behavior of -i/aAf^B-r is 
consistent with Tp in Fig. [2ja) in that the energy scale 
becomes smaller toward J C 2 , and has the asymmetry be- 
tween itinerant and localized sides. 

The magnetic response near the ILT shows a peculiar 
enhancement as shown in Fig. [2jd). We derive the local 
magnetic susceptibility defined by 



Xlo, 



= lim<#> 

fli-K) Hi 



(4) 



where Hi is the local magnetic field applied only at site 
i. Note that xi oc is independent of the site index even 
in the AFM phase. The results are shown in Fig. [2][d). 
For Jk < 0.15 and Jk > 0.25, xioc is almost the same 
for T = 0.0025 and T = 0.0010 shown in the figure. 
The smallness of xioc in the small Jk (< 0.15) region 
is due to the large magnetic moment, while in the large 
Jk (> 0.25) region, the smallness is due to the formation 
of Kondo singlets. 

Near the ILT point, the local magnetic susceptibility 
Xioc continuously increases with decreasing temperature, 
even though the magnetic moment is already large at 
J C 2 as shown in Fig. [2][b). The large fluctuation in the 
itinerant regime rapidly decreases on entering the local- 
ized side. We have checked (not shown) that uniform 
and staggered susceptibilities also show the rapid change 
at J C 2. The asymmetry between itinerant and localized 
sides is similar to Tp as shown in Fig.[2Ja), 

This anomalous magnetic fluctuation deep inside the 
magnetic phase originates from the formation of ex- 
tremely narrow electronic bands caused by the Kondo 
effect. The distinction between the AFM phases above 
and below J C 2 is most clearly visible in the single-particle 
spectral function defined by 



1 f°° dt 



-iEt 



(5) 



where 0(t) = e mt Oe- mt and f(E) = l/(e E / T + 1). 
The k dependence of A„(ek, E) in the DMFT enters only 
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FIG. 3: (Color online) Single-particle spectrum A a (sk, E) at 
T = 0.0025 for (a) J K = 0.125 and (b) J K = 0.225. The insets 
in (a) and (b) show the schematic illustrations for spectra in 
the localized (AFM) and itinerant (Para) limits, respectively. 



through £fc since the self energy has no wave-vector de- 
pendence. In the present notation, the Brillouin zone 
center and corner are located at Ek = —1/2 and Sk = 1/2, 
respectively. The energy £fe = gives the center of band 
folding by the cell doubling. Namely, £& = gives the 
boundary of the new Brillouin zone (BZ). Note that the 
spectrum (|5|) is independent of spin in the AFM phase. 

Figure [3] shows the single-particle spectrum at T = 
0.0025. In the localized regime shown in Fig. EJa), the 
Fermi surface is located near the new BZ boundary. The 
location of the minimum energy gap is the same as that in 
the localized limit Jk — > 0, as illustrated in the inset. In 
the itinerant regime, on the other hand, the Fermi surface 
is almost the same as in the paramagnetic heavy Fermi 
liquid as shown in Fig.^b) and its inset. Hence this state 
is classified as itinerant. In the itinerant regime, a very 
flat band appears at the Fermi level, while in the localized 
side the band is more dispersive at E = as shown in 
Fig. EJa) . This is related to the difference between the 
main components of the Fermi surface: mainly c electrons 
in the localized regime, but mainly / electrons in the 
itinerant side. It is this nearly flat band that causes the 
large magnetic fluctuation as seen in Fig. Efd). We note 
that the present model has no charge degrees of freedom 
for / electrons. Hence it is appropriate to interpret the 
emergent flat band in terms of coherent interaction of 
Kondo peaks, rather than direct cf hybridization. 

As shown in Fig. [IJa), a partly hybridized band struc- 
ture is observed even for Jk < J C 2, which belongs to our 
localized regime. This situation makes it necessary to 
define precisely the meaning of localized electrons. We 
rely on the character of the Fermi surface, since this dis- 
tinction is most practical and accessible by experiment. 
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FIG. 4: (Color online) Phase diagram of the Kondo- 

Heisenberg lattice in the Jh-Jk plane at T = 0.001. The 
abbreviations L- and I- mean localized and itinerant, respec- 
tively. At the critical value Jk = J C 2 for the ILT, the corre- 
sponding Jh mostly scales with the Kondo temperature Tk. 



Indeed, the ILT separates the two regimes with rather 
different physical properties as seen in Figs. [1] and [2] 

We have performed similar calculations for various val- 
ues of Jk and Jh, and obtained the phase diagram of the 
Kondo- Hcisenberg lattice. Figure @] shows the result at 
T = 0.001, which we regard as close to the ground state. 
Note, however, that we have not considered superconduc- 
tivity since the DMFT cannot deal with pairings other 
than the s-wave. The label L or I in the figure means the 
localized or itinerant magnetism, respectively. The blank 
circles and squares show the points where a discontinuous 
change takes place. At Jh = 0, which corresponds to the 
ordinary Kondo lattice, the AFM phase has an itinerant 
ground state regardless of the value of Jk (< J c i)- For 
small Kondo interaction, however, the itinerant AFM is 
fragile and easily turns into the localized one beyond a 
critical Jh (> 0). The critical value of Jh for the ILT 
with given Jk scales well with Tk as shown in Fig. [4] 
This clearly shows that the ILT in the AFM phase oc- 
curs if the Kondo temperature exceeds the Hcisenberg 
interaction. 

Next let us discuss the FM state for Jh < 0. For 
small enough | Jh | , the localized fcrromagnctism (L-FM) 
changes into itinerant AFM with increasing Jk , which is 
indicated by J C 4 in Fig. @] For larger Heisenberg inter- 
action with | Jh| > 0.07, the transition occurs between I- 
and L-FMs. There is another phase transition between 
I-FM and paramagnetic phases, and the corresponding 
value is indicated by J c3 in Fig. [4] The I-FM phase 
seems to persist up to large Kondo couplings, and can 
be understood as the spin selective Kondo insulator [22j . 
This state is interpreted as a mixture of fully polarized 
localized spins and local Kondo singlets. Here c electrons 
with the minority spin are insulating, while those with 



the majority spin are metallic. We note that the I- AFM 
cannot be understood in a manner similar to the I-FM, 
since the I- AFM phase is not connected to Jk — >• oo. 

Unlike the AFM, the itinerant and localized FMs can 
be connected either discontinuously or continuously. In 
the present case with n = 0.95, the L-FM changes to the 
I-FM by a first-order transition. As n decreases, how- 
ever, the first-order transition changes into a continu- 
ous crossover. We have actually observed a continuous 
change for n < 0.5. In Fig. [4j we have defined the L-FM 
state as the phase which is continuously connected to the 
localized limit with Jk = 0. 

Let us discuss relevance of our results to magnetic 
properties in real systems at low temperatures. With 
increasing pressure, the Kondo interaction will become 
larger in Ce and U systems. The arrows in Fig. |4] show 
this trend for Jk- On the other hand, it is not so simple 
to simulate the pressure dependence of Jh- Hence angle 
of the arrows do not have a microscopic basis. Let us first 
begin with the arrow (I) that simulates the behavior in 
CeRh x Coi_ x In5 [3[. Namely, decrease of x corresponds 
to the chemical pressure, and we expect the L-AFM to 
change into I-AFM and then to become paramagnetic. 
The magnetic structure in the I-AFM seems different 
from that in the L-AFM as observed in this compound 
[231 ] . With fine tuning of parameters, J c \ can be made 
very close to J C 2- Then the ILT occurs very close to the 

CP, which has been seen in CeRhIn,5 under pressure 
ll 
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Next we turn to the arrow (II) in Fig. |4] that simulates 
CcRu2(Si a: Gei_ 2 ;)2 |5J|. CcRu2Gc2 with x — is a local- 
ized ferromagnet, which changes into itinerant antiferro- 
magnet with increasing Si content, and finally becomes 
paramagnetic at x = 1. If Jh is more ferromagnetic, 
we obtain the arrow (III) in Fig. 0] that simulates UGe2 
under pressure, where the Fermi surface changes as the 
localized fcrromagnctism changes into the itinerant one 
inside the magnetic phase [6[. 

Another interesting example is found in Kondo insula- 
tors with magnetic order. Namely, CeRu2Alio is a heavy- 
electron antiferromagnet [25J, |26| , which is regarded as lo- 
cated close to the I-AFM and L-AFM boundary as illus- 
trated in Fig. |4] Slight substitution of Ru by Rh changes 
its itinerant character into a localized one [271 ] . Instead of 
Ru, we obtain CeOs2Al i0 [25|, which has a larger Kondo 
interaction, and should be located near J cl . On the other 
hand, CeFe2Alio [25| has the strongest cf interaction 
among the CeT2Alio family with T = Ru, Os, Fe and 
remains paramagnetic down to lowest accessible temper- 
atures. The crosses in Fig. [4] indicate locations of these 
compounds. Note that the values of Jh are uncertain 
and not important for sufficiently positive Jh- 

Finally, we discuss implication of the strong local mag- 
netic fluctuations near the ILT as shown in Fig. EJd). 
The entropy associated with the fluctuations should be 
removed at sufficiently low temperatures. It is tempting 



to interpret the superconductivity in CeRh^Coi-^Ins [3| 
as a result of releasing the entropy. In the boundary of L- 
FM and I-FM, the local fluctuation becomes also strong 
especially in the itinerant side. Hence the superconduc- 
tivity observed in UGe2 28j seems understandable by the 
same mechanism. 

In conclusion, we have identified the ILT inside ordered 
states of the Kondo-Heisenberg lattice by means of the 
DMFT. An extremely small energy scale appears near the 
transition point. It has been proposed that the observed 
phase diagrams in CeRh a: Coi_. r In5, CcRu^SizGei-a;^, 
UGe2 and CeT 2 Ali can be understood by the present 
model in a unified way. 
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